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Library Routine M 26 - 26k 



TITLE: 



TYPE: 

NUMBER OF WORDS: 

TEMPORARY STORAGE: 

DURATION: 

PARAMETERS: 



Eigenvalues and Eigenvectors of a Symmetric Jfatrix 

(SADOI Only) 

Complete program 

176 + OR 1) + (N 12) + (P 16) = 280 

0-30 

Approximately 5n^ milliseconds per iteration 

Four Parameters must be read in before the main routine 

and the matrix are read in, 

S3: OOF OOgF where 1 < g < 12 is the number of decimal 

digits to be printed. 

where n « order of matrix. 

where m = 2n if eigenvalues and eigenvectors 

both are desired; otherwise m - n. 
S6: OOF 00(311 + n(n+l)/2)F 
SIMBOUC ADDRESSES USED: (M26); (l); (2)$ (3); (V). 
LIBRARY ROUTINES USED: (Pl6); (N12); (Rl). 



3kz OOF OOnF 
S5: OOF OOmF 



LIMITATIONS: 



SCALING: 



If eigenvalues only are desired then n, the order of the 
matrix, must be < 37 • If eigenvectors are desired as 
well then n < 21. 
The elements of the matrix should be scaled so that 
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METHOD OF USE: 



The eigenvalues and eigenvectors will be scaled by the 
same factor. However, in order to prevent overflow the 
eigenvectors are actually scaled down by an additional 
factor of 10" 1 . 

1. Parameters must be read in first. 

2. Read in Program tape; machine will stop on a 2^(l0j) 

3. Insert matrix tape in reader, the elements of the 
matrix being punched as signed decimal fractions, i.e., 
as + or - followed by up to 12 decimal digits. As only 
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OUTPUT FJRMAT: 



half the off diagonal elements are needed (because 

of symmetry), the user need punch only a ; a ■■•-.' 

11 21* 22* 
*31' a 32' a 33* etc - 

i.e., we punch only /a. 

a 21 a 22 
a 31 a 32 a 33 





Following ajm an H terminating symbol must be punched. 

Start with black switch. 

If there has been an arithmetic failure then tvo F»s 

will be punched out and the machine will come to an 

OF stop at location (OLK)^. If the user, nevertheless, 

wishes to continue the computations he may do so by 

raising and lowering the white switch. 

If there has been no arithmetic failure then the first 

thing punched out will be an integer - the number of 

iterations required for convergence. This will be 

followed by the column of eigenvalues. The machine 

will then come to a 3* stop at (121)^. If eigenvectors 

are also desired then a black switch start will cause 

the eigenvectors to be punched as n columns of n components 

each (the eigenvectors being separated from one another 

by an extra line feed and carriage return). The machine 

will then come to an OF stop at (132) 

SD 

For a complete description of the mathematical method 

see the write-up of M?5. 
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